Linkage disequilibrium assessment via log-linear modeling of SNP haplotype frequencies

Analyses of high-density single-nucleotide polymorphism (SNP) data, such as genetic mapping and linkage disequilibrium (LD) studies, require phase-known haplotypes to allow for the correlation between tightly linked loci. However, current SNP genotyping technology cannot determine phase, which must be inferred statistically. In this paper, we present a new Bayesian Markov chain Monte Carlo (MCMC) algorithm for population haplotype frequency estimation, particulary in the context of LD assessment. The novel feature of the method is the incorporation of a log-linear prior model for population haplotype frequencies. We present simulations to suggest that 1) the log-linear prior model is more appropriate than the standard coalescent process in the presence of recombination (>0.02cM between adjacent loci), and 2) there is substantial inflation in measures of LD obtained by a "two-stage" approach to the analysis by treating the "best" haplotype configuration as correct, without regard to uncertainty in the recombination process. Genet Epidemiol 25:106-114, 2003. (C) 2003 Wiley-Liss, Inc.

Development of a dense SNP-based linkage map of an apple rootstock progeny using the Malus Infinium whole genome genotyping array

Background A whole-genome genotyping array has previously been developed for Malus using SNP data from 28 Malus genotypes. This array offers the prospect of high throughput genotyping and linkage map development for any given Malus progeny. To test the applicability of the array for mapping in diverse Malus genotypes, we applied the array to the construction of a SNPbased linkage map of an apple rootstock progeny. Results Of the 7,867 Malus SNP markers on the array, 1,823 (23.2 %) were heterozygous in one of the two parents of the progeny, 1,007 (12.8 %) were heterozygous in both parental genotypes, whilst just 2.8 % of the 921 Pyrus SNPs were heterozygous. A linkage map spanning 1,282.2 cM was produced comprising 2,272 SNP markers, 306 SSR markers and the S-locus. The length of the M432 linkage map was increased by 52.7 cM with the addition of the SNP markers, whilst marker density increased from 3.8 cM/marker to 0.5 cM/marker. Just three regions in excess of 10 cM remain where no markers were mapped. We compared the positions of the mapped SNP markers on the M432 map with their predicted positions on the ‘Golden Delicious’ genome sequence. A total of 311 markers (13.7 % of all mapped markers) mapped to positions that conflicted with their predicted positions on the ‘Golden Delicious’ pseudo-chromosomes, indicating the presence of paralogous genomic regions or misassignments of genome sequence contigs during the assembly and anchoring of the genome sequence. Conclusions We incorporated data for the 2,272 SNP markers onto the map of the M432 progeny and have presented the most complete and saturated map of the full 17 linkage groups of M. pumila to date. The data were generated rapidly in a high-throughput semi-automated pipeline, permitting significant savings in time and cost over linkage map construction using microsatellites. The application of the array will permit linkage maps to be developed for QTL analyses in a cost-effective manner, and the identification of SNPs that have been assigned erroneous positions on the ‘Golden Delicious’ reference sequence will assist in the continued improvement of the genome sequence assembly for that variety.

Evaluation of the use of high-density SNP genotyping to implement UPOV Model 2 for DUS testing in barley

Developments in high-throughput genotyping provide an opportunity to explore the application of marker technology in distinctness, uniformity and stability (DUS) testing of new varieties. We have used a large set of molecular markers to assess the feasibility of a UPOV Model 2 approach: “Calibration of threshold levels for molecular characteristics against the minimum distance in traditional characteristics”. We have examined 431 winter and spring barley varieties, with data from UK DUS trials comprising 28 characteristics, together with genotype data from 3072 SNP markers. Inter varietal distances were calculated and we found higher correlations between molecular and morphological distances than have been previously reported. When varieties were grouped by kinship, phenotypic and genotypic distances of these groups correlated well. We estimated the minimum marker numbers required and showed there was a ceiling after which the correlations do not improve. To investigate the possibility of breaking through this ceiling, we attempted genomic prediction of phenotypes from genotypes and higher correlations were achieved. We tested distinctness decisions made using either morphological or genotypic distances and found poor correspondence between each method.

Analysis of DNA polymorphism in ancient barley herbarium material: validation of the KASP SNP genotyping platform

The availability of crop specimens archived in herbaria and old seed collections represent valuable resources for the analysis of plant genetic diversity and crop domestication. The ability to extract ancient DNA (aDNA) from such samples has recently allowed molecular genetic investigations to be undertaken in ancient materials. While analyses of aDNA initially focused on the use of markers which occur in multiple copies such as the internal transcribed spacer region (ITS) within ribosomal DNA and those requiring amplification of short DNA regions of variable length such as simple sequence repeats (SSRs), emphasis is now moving towards the genotyping of single nucleotide polymorphisms (SNPs), traditionally undertaken in aDNA by Sanger sequencing. Here, using a panel of barley aDNA samples previously surveyed by Sanger sequencing for putative causative SNPs within the flowering-time gene PPD-H1, we assess the utility of the Kompetitive Allele Specific PCR (KASP) genotyping platform for aDNA analysis. We find KASP to out-perform Sanger sequencing in the genotyping of aDNA samples (78% versus 61% success, respectively), as well as being robust to contamination. The small template size (≥46 bp) and one-step, closed-tube amplification/genotyping process make this platform ideally suited to the genotypic analysis of aDNA, a process which is often hampered by template DNA degradation and sample cross-contamination. Such attributes, as well as its flexibility of use and relatively low cost, make KASP particularly relevant to the genetic analysis of aDNA samples. Furthermore, KASP provides a common platform for the genotyping and analysis of corresponding SNPs in ancient, landrace and modern plant materials. The extended haplotype analysis of PPD-H1 undertaken here (allelic variation at which is thought to be important for the spread of domestication and local adaptation) provides further resolution to the previously identified geographic cline of flowering-time allele distribution, illustrating how KASP can be used to aid genetic analyses of aDNA from plant species. We further demonstrate the utility of KASP by genotyping ten additional genetic markers diagnostic for morphological traits in barley, shedding light on the phenotypic traits, alleles and allele combinations present in these unviable ancient specimens, as well as their geographic distributions.

A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.)

Faba bean (Vicia faba L.) is a globally important nitrogen-fixing legume, which is widely grown in a diverse range of environments. In this work, we mine and validate a set of 845 SNPs from the aligned transcriptomes of two contrasting inbred lines. Each V. faba SNP is assigned by BLAST analysis to a single Medicago orthologue. This set of syntenically anchored polymorphisms were then validated as individual KASP assays, classified according to their informativeness and performance on a panel of 37 inbred lines, and the best performing 757 markers used to genotype six mapping populations. The six resulting linkage maps were merged into a single consensus map on which 687 SNPs were placed on six linkage groups, each presumed to correspond to one of the six V. faba chromosomes. This sequence-based consensus map was used to explore synteny with the most closely-related crop species, lentil, and the most closely related fully sequenced genome, Medicago. Large tracts of uninterrupted colinearity were found between faba bean and Medicago, making it relatively straightforward to predict gene content and order in mapped genetic interval. As a demonstration of this, we mapped a flower colour gene to a 2 cM interval of Vf chromosome 2 which was highly collinear with Mt3. The obvious candidate gene from 77 gene models in the collinear Medicago chromosome segment was the previously characterized MtWD40-1 gene (Mt3g092830, Mt3g092840) controlling anthocyanin production in Medicago and re-sequencing of the Vf orthologue showed a putative causative deletion of the entire 5’ end of the gene.

Flanking SNP markers for vicine–convicine concentration in faba bean (Vicia faba L).

The pyrimidine glycosides, vicine and convicine, limit the use of faba bean (Vicia faba L.) as food and feed. A single recessive gene, vc-, is responsible for a lowered vicine–convicine concentration. The biosynthetic pathway of these closely related compounds is not known, and the nearest available markers are several cM away from vc-. Improved markers would assist breeding and help to identify candidate genes. A segregating population of 210 F5 recombinant inbred lines was developed from the cross of Mélodie/2 (low vicine–convicine) × ILB 938/2 (normal vicine–convicine), and vicine–convicine concentrations were determined twice on each line. The population was genotyped with a set of 188 SNPs. A strong, single QTL for vicine–convicine concentration was identified on chromosome I, flanked by markers 1.0 cM away on one side and 2.6 cM on the other. The interval defined by these markers in the model species Medicago truncatula includes about 340 genes, but no candidate genes were identified. Further fine mapping should lead to the identification of tightly linked markers as well as narrowing down the search for candidate regulatory or biosynthetic genes which could underlie the vc- locus.

Predicting fatty acid profiles in blood based on food intake and the FADS1 rs174546 SNP

SCOPE: A high intake of n-3 PUFA provides health benefits via changes in the n-6/n-3 ratio in blood. In addition to such dietary PUFAs, variants in the fatty acid desaturase 1 (FADS1) gene are also associated with altered PUFA profiles. METHODS AND RESULTS: We used mathematical modelling to predict levels of PUFA in whole blood, based on MHT and bolasso selected food items, anthropometric and lifestyle factors, and the rs174546 genotypes in FADS1 from 1,607 participants (Food4Me Study). The models were developed using data from the first reported time point (training set) and their predictive power was evaluated using data from the last reported time point (test set). Amongst other food items, fish, pizza, chicken and cereals were identified as being associated with the PUFA profiles. Using these food items and the rs174546 genotypes as predictors, models explained 26% to 43% of the variability in PUFA concentrations in the training set and 22% to 33% in the test set. CONCLUSIONS: Selecting food items using MHT is a valuable contribution to determine predictors, as our models' predictive power is higher compared to analogue studies. As unique feature, we additionally confirmed our models' power based on a test set.

p22phox C242T SNP inhibits inflammatory oxidative damage to endothelial cells and vessels

Background— T NADPH oxidase, by generating reactive oxygen species, is involved in the pathophysiology of many cardiovascular diseases and represents a therapeutic target for the development of novel drugs. A single-nucleotide polymorphism (SNP) C242T of the p22phox subunit of NADPH oxidase has been reported to be negatively associated with coronary heart disease (CHD) and may predict disease prevalence. However, the underlying mechanisms remain unknown. Methods and Results— Using computer molecular modelling we discovered that C242T SNP causes significant structural changes in the extracellular loop of p22phox and reduces its interaction stability with Nox2 subunit. Gene transfection of human pulmonary microvascular endothelial cells showed that C242T p22phox reduced significantly Nox2 expression but had no significant effect on basal endothelial O2.- production or the expression of Nox1 and Nox4. When cells were stimulated with TNFα (or high glucose), C242T p22phox inhibited significantly TNFα-induced Nox2 maturation, O2.- production, MAPK and NFκB activation and inflammation (all p<0.05). These C242T effects were further confirmed using p22phox shRNA engineered HeLa cells and Nox2-/- coronary microvascular endothelial cells. Clinical significance was investigated using saphenous vein segments from non CHD subjects after phlebectomies. TT (C242T) allele was common (prevalence of ~22%) and compared to CC, veins bearing TT allele had significantly lower levels of Nox2 expression and O2.- generation in response to high glucose challenge. Conclusions— C242T SNP causes p22phox structural changes that inhibit endothelial Nox2 activation and oxidative response to TNFα or high glucose stimulation. C242T SNP may represent a natural protective mechanism against inflammatory cardiovascular diseases.

Family-based association analysis with ordered categorical phenotypes, covariates and interactions

Genetic association analyses of family-based studies with ordered categorical phenotypes are often conducted using methods either for quantitative or for binary traits, which can lead to suboptimal analyses. Here we present an alternative likelihood-based method of analysis for single nucleotide polymorphism (SNP) genotypes and ordered categorical phenotypes in nuclear families of any size. Our approach, which extends our previous work for binary phenotypes, permits straightforward inclusion of covariate, gene-gene and gene-covariate interaction terms in the likelihood, incorporates a simple model for ascertainment and allows for family-specific effects in the hypothesis test. Additionally, our method produces interpretable parameter estimates and valid confidence intervals. We assess the proposed method using simulated data, and apply it to a polymorphism in the c-reactive protein (CRP) gene typed in families collected to investigate human systemic lupus erythematosus. By including sex interactions in the analysis, we show that the polymorphism is associated with anti-nuclear autoantibody (ANA) production in females, while there appears to be no effect in males.

A functional genomics approach reveals novel quantitative trait loci associated with platelet signalling pathways

Platelet response to activation varies widely between individuals but shows interindividual consistency and strong heritability. The genetic basis of this variation has not been properly explored. We therefore systematically measured the effect on function of sequence variation in 97 candidate genes in the collagen and adenosine-diphosphate (ADP) signaling pathways. Resequencing of the genes in 48 European DNA samples nearly doubled the number of known single nucleotide polymorphisms (SNPs) and informed the selection of 1327 SNPs for genotyping in 500 healthy Northern European subjects with known platelet responses to collagen-related peptide (CRP-XL) and ADP. This identified 17 novel associations with platelet function (P < .005) accounting for approximately 46% of the variation in response. Further investigations with platelets of known genotype explored the mechanisms behind some of the associations. SNPs in PEAR1 associated with increased platelet response to CRP-XL and increased PEAR1 protein expression after platelet degranulation. The minor allele of a 3' untranslated region (UTR) SNP (rs2769668) in VAV3 was associated with higher protein expression (P = .03) and increased P-selectin exposure after ADP activation (P = .004). Furthermore the minor allele of the intronic SNP rs17786144 in ITPR1 modified Ca2+ levels after activation with ADP (P < .004). These data provide novel insights into key hubs within platelet signaling networks.

A novel variant on chromosome 7q22.3 associated with mean platelet volume, counts, and function

Mean platelet volume (MPV) and platelet count (PLT) are highly heritable and tightly regulated traits. We performed a genome-wide association study for MPV and identified one SNP, rs342293, as having highly significant and reproducible association with MPV (per-G allele effect 0.016 +/- 0.001 log fL; P < 1.08 x 10(-24)) and PLT (per-G effect -4.55 +/- 0.80 10(9)/L; P < 7.19 x 10(-8)) in 8586 healthy subjects. Whole-genome expression analysis in the 1-MB region showed a significant association with platelet transcript levels for PIK3CG (n = 35; P = .047). The G allele at rs342293 was also associated with decreased binding of annexin V to platelets activated with collagen-related peptide (n = 84; P = .003). The region 7q22.3 identifies the first QTL influencing platelet volume, counts, and function in healthy subjects. Notably, the association signal maps to a chromosome region implicated in myeloid malignancies, indicating this site as an important regulatory site for hematopoiesis. The identification of loci regulating MPV by this and other studies will increase our insight in the processes of megakaryopoiesis and proplatelet formation, and it may aid the identification of genes that are somatically mutated in essential thrombocytosis. (Blood. 2009; 113: 3831-3837)

Sequentially testing for a gene-drug interaction in a genomewide analysis

Assaying a large number of genetic markers from patients in clinical trials is now possible in order to tailor drugs with respect to efficacy. The statistical methodology for analysing such massive data sets is challenging. The most popular type of statistical analysis is to use a univariate test for each genetic marker, once all the data from a clinical study have been collected. This paper presents a sequential method for conducting an omnibus test for detecting gene-drug interactions across the genome, thus allowing informed decisions at the earliest opportunity and overcoming the multiple testing problems from conducting many univariate tests. We first propose an omnibus test for a fixed sample size. This test is based on combining F-statistics that test for an interaction between treatment and the individual single nucleotide polymorphism (SNP). As SNPs tend to be correlated, we use permutations to calculate a global p-value. We extend our omnibus test to the sequential case. In order to control the type I error rate, we propose a sequential method that uses permutations to obtain the stopping boundaries. The results of a simulation study show that the sequential permutation method is more powerful than alternative sequential methods that control the type I error rate, such as the inverse-normal method. The proposed method is flexible as we do not need to assume a mode of inheritance and can also adjust for confounding factors. An application to real clinical data illustrates that the method is computationally feasible for a large number of SNPs. Copyright (c) 2007 John Wiley & Sons, Ltd.

Influence of apoA-V gene variants on postprandial triglyceride metabolism: impact of gender

Although apolipoprotein AN (apoA-V) polymorphisms have been consistently associated with fasting triglyceride (TG) levels, their impact on postprandial lipemia remains relatively unknown. In this study, we investigate the impact of two common apoA-V polymorphisms (-1131 T>C and S19W) and apoA-V haplotypes on fasting and postprandial lipid metabolism in adults in the United Kingdom (n = 259). Compared with the wild-type TT, apoA-V -1131 TC heterozygotes had 15% (P = 0.057) and 21% (P = 0.002) higher fasting TG and postprandial TG area under the curve (AUC), respectively. Significant (P = 0.038) and nearly significant (P = 0.057) gender X genotype interactions were observed for fasting TG and TG AUC, with a greater impact of genotype in males. Lower HDL-cholesterol was associated with the rare TC genotype (P = 0.047). Significant linkage disequilibrium was found between the apoA-V -1131 T>C and the apoC-III 3238 C>G variants, with univariate analysis indicating an impact of this apoC-III single nucleotide polymorphism (SNP) on TG AUC (P = 0.015). However, in linear regression analysis, a significant independent association with TG AUC (P = 0.007) was only evident for the apoA-V -1131 T>C SNP, indicating a greater relative importance of the apoA-V genotype.

The impact of age on the postprandial vascular response to a fish oil-enriched meal

Although chronic fish oil intervention had been shown to have a positive impact on vascular reactivity, very little is known about their acute effects during the postprandial phase. Our aim was to examine the impact of a fish oil-enriched test meal on postprandial vascular reactivity in healthy younger ( < 50 years) v. older ( ≥ 50 years) men. Vascular reactivity was measured at baseline (0 h), 2 and 4 h after the meal by laser Doppler iontophoresis and blood samples taken at 0 and 4 h for the measurement of plasma lipids, total nitrite, glucose and insulin. Acetylcholine- (ACh, endothelial-dependent vasodilator) and sodium nitroprusside (SNP, endothelial-independent vasodilator)-induced reactivities were greater at 4 h than at baseline or 2 h in the younger men (P < 0·04). These changes were not observed in the older men. Comparison of the male groups revealed significantly greater responses to ACh (P = 0·006) and SNP (P = 0·05) at 4 h in the younger compared with the older males. Postprandial NEFA concentrations were also greater at 4 h in the younger compared with the older men (P = 0·005), with no differences observed for any of the other analytes. Multiple regression analysis revealed age to be the most significant predictor of both ACh and SNP induced reactivity 4 h after the meal. In conclusion, the ingestion of a meal enriched in fish oil fatty acids was shown to improve postprandial vascular reactivity at 4 h in our younger men, with little benefit evident in our older men.

Gene-nutrient interactions in the metabolic syndrome: single nucleotide polymorphisms in ADIPOQ and ADIPOR1 interact with plasma saturated fatty acids to modulate insulin resistance

Background: Progression of the metabolic syndrome (MetS) is determined by genetic and environmental factors. Gene-environment interactions may be important in modulating the susceptibility to the development of MetS traits. Objective: Gene-nutrient interactions were examined in MetS subjects to determine interactions between single nucleotide polymorphisms (SNPs) in the adiponectin gene (ADIPOQ) and its receptors (ADIPOR1 and ADIPOR2) and plasma fatty acid composition and their effects on MetS characteristics. Design: Plasma fatty acid composition, insulin sensitivity, plasma adiponectin and lipid concentrations, and ADIPOQ, ADIPOR1, and ADIPOR2 SNP genotypes were determined in a cross-sectional analysis of 451 subjects with the MetS who participated in the LIPGENE (Diet, Genomics, and the Metabolic Syndrome: an Integrated Nutrition, Agro-food, Social, and Economic Analysis) dietary intervention study and were repeated in 1754 subjects from the LIPGENE-SU.VI.MAX (SUpplementation en VItamines et Mineraux AntioXydants) case-control study (http://www.ucd.ie/lipgene). Results: Single SNP effects were detected in the cohort. Triacylglycerols, nonesterified fatty acids, and waist circumference were significantly different between genotypes for 2 SNPs (rs266729 in ADIPOQ and rs10920533 in ADIPOR1). Minor allele homozygotes for both of these SNPs were identified as having degrees of insulin resistance, as measured by the homeostasis model assessment of insulin resistance, that were highly responsive to differences in plasma saturated fatty acids (SFAs). The SFA-dependent association between ADIPOR1 rs10920533 and insulin resistance was replicated in cases with MetS from a separate independent study, which was an association not present in controls. Conclusions: A reduction in plasma SFAs could be expected to lower insulin resistance in MetS subjects who are minor allele carriers of rs266729 in ADIPOQ and rs10920533 in ADIPOR1. Personalized dietary advice to decrease SFA consumption in these individuals may be recommended as a possible therapeutic measure to improve insulin sensitivity. This trial was registered at clinicaltrials.